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ABSTRACT 

We show that time variations in the UV ionizing continuum of quasars, on 
scales of ~1 year, affect the dynamic structure of the plasmas responsible for low 
ionization broad absorption lines. Variations of the ionizing continuum produce 
non-equilibrium photoionization conditions over a significant fraction of the ab- 
sorbing clouds and supersonically moving ionization fronts. When the flux drops 
the contraction of the ionized region drives a supersonic cooling front towards the 
radiation source and a rarefaction wave in the opposite direction. The pressure 
imbalance is compensated by an increased speed of the cool gas relative to the 
front. When the flux recovers the cool gas is re- ionized and re- heated by a su- 
personic ionization front traveling away from the radiation source and a forward 
shock is created. The reheated clouds equilibrate to a temperature of ~ 10"^ K 
and are observed to have different radial velocities than the main cloud. Such 
fragmentation seems consistent with the multicomponent structure of troughs 
seen in some objects. The velocity differences measured among various compo- 
nents in the quasars QSO 2359-1241 and SDSS J0318-0600 can be reproduced 
by our model if strong magnetic fields (~10 mG) are present within the clouds. 



1. Introduction 



Gas outflows a potential source of feedback that regulates the evolution of the central 
black hole, the host galaxy, and the inter galactic medium (e.g. Elvis 2006). In the rest-frame 



UV spectra of 10 - 20% of all quasars (IKnigge et al.l l2008l : iFoltz et al.l 119871 ). outflows are 



manifested in blueshifted Broad Absorption Lines (BAL). The most common of these are 



Wevmann 


1995; 


Turnshek 


1995) 



Unfortunately, their saturated nature limits their utility in determining the physical charac- 
teristics of the outflow. A subset of BAL quasars, known as LoBALs, also show complexes 
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of narrower (on the order of ~ 100 km s~^) absorption features from low ionization species 
such as Mg II, Al II, and Si II. Among these, objects that include troughs of metastable Fe II 
or Fe III are specifically called FeLoBALs. These are useful because the column densities 
of metastable levels offer valuable diagnostics of the outflow (Gabel et al. 2006; Arav et 
al. 2005; and references therein) . These studies place the gas respons ible for the FeLoBAL 
at distances between ~ 1 kpc (IKorista et al.l l2008l : iMoe et al.l l2009l ) and up to ~10 kpc 



flDunn et al.ll2010f l. 



We studied the various kinematic components in the FeLoBAL troughs of two quasars 
QSO 2359-1241 and SDSS J0318-0600, z =0.868 and 1.967 respectively feautista et aL 
20101 ). It was found that among the various kinematic components the highest column 



density component of each system also had the highest particle density, while all other 
components have densities of about one quarter of that of the main component. Moreover, all 
the components seem to be located at the same distance from the central source, suggesting 



that t hey are all related to each other. This is also consistent with the finding of IVoit et al. 



(119931 ) that BALs and LoBALs could be explained simultaneously by ablation of either a 
single large cloud or many scattered clouds. 

Another characteristic of FeLoBALs is that, within the sample of known objects, they 
are borderline density-bounded. FeLoBAL quasars are seen in absorption from species that 
arise from the ionization front (IF) of the cloud (e.g. Si II and Fe II) while atoms from a 
neutral cold region, e.g. C I or Fe I, are not seen. This seems to point out some relationship 
between the ionization front and the nature of the clouds. The location of the IF of a cloud 
of a given density depends on the flux of hydrogen ionizing photons, particularly in the ~1 
to 4 Ryd energy range, from the radiation source. However, quasar con tinua are known t o 
vary in the optical and UV within time scales of the order of one year (IKaspi et al.l 120001 ). 
Hence, it is important to understand the effect of the varying radiation field on the structure 
of the FeLoBAL. 

In this letter we show that when the flux that ionizes a FeLoBAL cloud drops the 
depth of the ionized region often shrinks at supersonic speeds followed by a cooling front, 
as explained in Sections 2 and 3. In section 4 we show that the supersonic cooling front 
creates a rarefaction wave, traveling in the opposite direction to the front, which results in 
fragmentation of the cloud. Under some conditions, such a fragmentation can explain the 
observed structure of FeLoBALs. 
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2. Magnitude scales and the effect of flux variations on FeLoBALs 

Let us start with the simple picture illustrated in Fig. [TJ Here, r is the depth of the 
neutral region of the cloud and it is measured from the far end of the cloud towards the 
central source. For a given luminosity of the source and a cloud of constant particle density 
Uh moving in a medium of negligible density one can draw an imaginary radius within which 
the cloud is fully ionized. This is 

(sis)"'' 

where Qh is rate of hydrogen ionizing photons from the central source, a is the recombination 
rate coefficient, and / is the depth of the cloud. For a cloud located beyond Rif the ionized 
depth (/ — r) is related to the physical distance {R) as R^{1 — r) = Q „ / {Ann^jO) . Thus, 
FeLoBAL outflows that result from the formation of an IF, but without an extended neutral 
region, must be located right at Rif. Then, a drop AQ^ in the ionization flux, i.e. Q'j^ = 
Qh — ^Qh, will create a region of depth r = I x AQh/Q^ of either neutral or recombining 
plasma. Here, / ~ Aj^/n^, where is total column density of the cloud and has a typical 
value of ~ lO^'^'^ cm~^. In deriving this equation we assumed that / << Rif and r « I. 
In QSO 2359-1241 and SDSS J0318-0600 n„ are ~ 10^-^ cm'^ and ~ lO^-^ cm'^ and the 
sizes of the absorbers are 10^^-10^^ cm (i.e. ~0.01 to 0.1 light-years). The change in the 
ionization of the gas is not instantaneous, but delayed by the recombination time scale of the 
plasma = Ijn^a^ with a~ 2 x 10~^^ s~^cm~^ at Te ~ 10"^ K. Then, r„ 10 yrs in QSO 
2359-1241 and 100 yrs in SDSS J0318-0600. For example, in the case that Qh is abruptly 
shut off, the size of the ionized region is 

l^r{t) = le-'l^% (2) 

where t is the time measured from the moment that the flux is turned off. Thus, the depth 
of the ionized region shrinks with a speed 

d(l-r(t)) I 

Vif = ^ ^ — fort « r,. (3) 

at Ta 

Thus, from the derivation of / and above fj/ ~ Aj^a 600 km s~^, which is much faster 
than the speed of sound in the plasma. 

The width of quasi-static IFs is of the order of the mean free path for ionizing radiation, 
i.e. X/tih Oh, where s:^ 6 x 10 cm^ is the photoionization cross section. Thus, the IF's 
in QSO 2359-1241 and SDSS J03 18-0600 are on the order of lO^^-lO^^ cm, which are much 
smaller than the distance on which the IF retrats for flux variations on ~ 1 year time scale. 
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Still, these continuum variations much shorter than the recombination time will broaden 
the IF creating an extended recombining region. But, as the IF departs from ionization 
equilibrium it also departs from thermal equilibrium. In a static IF the temperature drops 
to a much lower value for a given value of the optical depth to heating photons. The optical 
depth is 

r 3 max 

r = cth n foods', (4) 



with rijjo the neutral hydrogen density. Let us assume, for instance, a static IF of the 
form n„o = Uh x (1 — e~^^^). As the ionized region shrinks the characteristic depth of 
the front increases to 6' > 6. At an early time t after an instantaneous retreat of the IF 
n^o t/ta- Thus, for an IF brodening by a factor u it will take a time t ^ t/v for the 
cooling front to start moving with the retreating IF. Then, if the value of r for which the 
temperature drops remains constant the depth within the IF at which this happens, Smaxi 
will be reduced with respect to 5' . For example, if the change in temperature occurred at 
Smax/S = 1 in the static IF, Smax/5' = 0.68 for 6'/ 6 = 2, Smax/S' = 0.41 for 6'/ 6 = 5, and 
Smax/S' = 0.28 for 6'/ 6 = 10. In the region of new lower equilibrium temperature cooling 
of the plasma proceeds muc h faster than r ecombination, because the cooling time scale is 



only tt ^ 2 X lO^/rin years (jSpitzer Ill998l ). In conclusion, the drop in temperature closely 
follows the supersonic shrinking of the ionized region despite the broadening IF. Note that 
a drop in temperature makes recombination speed up, because a increases as Te ^^"^ as the 
electron temperature decreases. 

What happens when the flux increases again? As the ionizing flux rises the neutral 
fragment of the cloud will be re- ionized. Thus, an IF will travel through the neutral fragment, 
in the direction away from the ionizing source. This process is governed by the equation 

4^ = (5) 

The speed of the IF in this case is proportional by the change in radiative flux and can be 
highly supersonic, in which case it will send shocks ahead of the IF and rarefaction waves in 
the opposite direction. 

In conclusion, time variations of the ionizing flux are expected to have two kinds of 
effects: (1) In the region neighboring the IF the ionization varies with time and departure 
from photoionization equilibrium could occur. This effect is out of the scope of this paper, 
but will be the subject of future investigation. (2) The varying ionizing flux will create 
supersonic cooling and heating fronts traveling through the ionized fraction of the cloud and 
the neutral region. These will be accompanied by rarefaction waves and shocks. This process 
is investigated in Section 4 of this letter, but first we will derive a more accurate description 
of the dynamics of the IF. 
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Fig. 1. — Illustration of a photoionized cloud of depth / traveling away from the light source. 
The upper cloud is right at Rif and the whole depth / is ionized. In the lower cloud only a 
depth / — r is ionized due to either a larger physical distance or a drop in Qh- 
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3. Dynamics of ionization fronts within outflowing clouds 



The IF within a cloud is defined as the region where the ionization drops to near zero. 
Thus, we solve for rte in the time dependent photoionization equation 



~dt 



+ V.{neu) = r - A, 



(6) 



Where t is time, u is the velocity of the gas, and F and A are the ionization and recombination 
rates per volume. 



-V.$ 



H 



V.- 



(7) 



A = a„nl 



where t{x) is the optical depth within the cloud up to the depth x. 

Integrating Eqn. 12] in the co-moving frame of the IF [u = 0) over the ionized volume 
V = VLR'^il — r), with Vt the covering solid angle, one finds a differential equation for r{t) 



d{l-r{t)) 
Jt 



/ - r{t) 



(9) 



For gas clouds much smaller than the distance to the ionizing source, i.e. / <C R, one can 
assume constant density during a crossing time of the cloud (~ l/v). Then, 



/ - r{t) = /e-*/- [ ^ f Quity^-dt' + 1 
jaQuKy^) Jo 



(10) 



To understand this solution it is useful to check a couple of cases: 

1) Qnit) varying as a step function. Let us assume that at t = the absorber just reached 
Rif when the continuum flux was Qo- Thus, the cloud is fully ionized and a drop in Qh will 
produce an IF. If the absorber were closer than Rif it could remain fully ionized even when 
Qh is reduced. If the flux drops to ^Qo at t = for t > one gets 

l-r{t) = lx [e + (l-Oe-*/"'^] (11) 

and 

vjF = {I ~ 0-e-'^^\ (12) 

which is consistent with the result of equation (4) and yields supersonic vjp for flux variations 
of ~ 5% or greater, i.e. ^ < 0.95. 
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Fig. 2. — Flux, r(t), and vipit) for a sample of PG quasars observed by lKaspi et al.l (l2000l ). 
The results are for PG0844 (left panels), PG1226 (middle panels), and PG1411 (right panels). 
The calculations assume Tq, = 120 yrs and / = 10^^ cm. The horizontal dotted lines in the 
viF VS. t plots delimit the approximate range for subsonic motion. 
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2) Quit) from observed quasar continua. Here we compute the io nization of the cloud vs. 
time for selected quasar continua observed by iKaspi et al.l (I2OOOI ) using Eqn.(lO). Fig. [2] 
shows the observer flux, the IF location with respect to the total depth of the cloud I and 
vjF for three quasars. These objects are chosen out of a sample of 17 quasars observed by 
Kaspi et al. because these illustrate the types of dynamical behaviors of the IF. Because 
the calculations are done for illustrative purposes only we disregard the uncertainties on 
the measured fluxes. Notice, for example, that in PG0844 the r{t)/l ratio is generally 
decreasing with time because the average ionizing continuum is lower than at the initial 
time of observation. This overall behavior is somewhat arbitrary as it depends on the initial 
observation. What is more meaningful, though, is the instantaneous velocity of the IF front, 
computed as Ar/At, and this is seen to reach highly supersonic values. This is because the 
flux levels in PG0844, PG1226, and PG1411 vary by ~40%, ~50%, and ~30% with respect 
to the value at t = in less than a year. Clearly, this determination of velocities is sensitive 
to the observational errors as well as the sampling of the quasar luminosity. Hence, we make 
no definitive claims on specific objects. However, the prediction that highly supersonic fronts 
are expected is supported in every object of the available sample. 



4. Supersonic ionization fronts and cooling fronts 

The supersonic IF is accompanied by a cooling front that leaves a pressure step behind. 
This creates a rearward- facing, i.e. in the opposite di r ection to the IF, rarefaction w ave and 
a forward-facing shock (jStrachan and Ahlborru Il975l : lAhlborn and Stracharu 119731 ). Under 
homogene ous conditions and for a n infinite flow stream these waves and shocks will form pe- 
riodically (jSchnerr fc Adamlll997l ). As discussed before, the fragments once created become 
re-ionized and re-heated when the flux recovers. Then one can study the jump conditions 
between the main cloud and the fragments assuming constant temperature. 

In one dimension the jump conditions across the rarefaction wave are 



PlWi = P2W2 



(13) 



Piwl + Pi + -i = p2wl + P2 + ^ 



(14) 



BiWi = B2W2 



(15) 



where pi is the mass density of the gas, w is the bulk velocity of the gas in the co-moving 
frame of the wave {w = v — Vwave), P is the thermal pressure, and B is the strength of the 
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magnetic field. Because the jump is essentially isothermal we adopt = Pi/ pi = P2/P2, 
and from Eqns. (13) - (15) one gets 



Pi 



(16) 



where uai = Bi/iSnpiY^"^ is the Alfven velocity. Thus, the difference of velocities between 
the fragments and the main cloud is 

1/2 



V2 — vi = W2 — wi ±20 X 



VAl 



4 \ 13 km s 



-1 



(17) 



where we have used typical conditions pxj P2 = 4 and a ^ 13 km s^^. This result shows that 
rarefaction waves, originated by supersonic IF, yield lower density fragments that detach 
from the main cloud with relative velocities of tens of km s~^ in absence of magnetic fields 
and much larger velocities when these are present. 

There is evidence suggesting the presence of magnetic fields in BALs. Large trough 
widths, of the order of 100 km/s, are observed in FeLoBALS, with the largest widths being 
generally associated with the main absorber (see Table 1 of Bautista et al. 2010). The 
widths are much greater than expected from thermal broadening and seem to indicate su- 
personic turbulent velocities which could only be sustained by equally fast Alfven waves. 
The alternative scenario that observed trough widths could be mostly due to bulk velocity 
gradients along the line of sight is inconsistent with the fact that clouds have been able to 
travel ~ 1 — 10 kpc away from the AGN without diffusing into the interstellar medium. 
Moreover, it is difficult to understand how bulk velocity gradients could yield similar trough 
profiles for ions with different ionization stages (but see Arav et al. 1999). 



Strong magnetic fields in BAL clouds have been proposed in the past. Ide Kool and Begelman 



( I1995I ) pointed out that magnetic pressures in BAL clouds must exceed thermal and radiative 
pressures to confine the clouds when accelerated up to thousands of km s~^. Their models 
for acceleration of BAL clou ds yield magnetic field strengths of the order of 10 mG. Further, 
Furlanetto and Loeb I (l200l[ ) argue that BAL quasar outflows could carry a large fraction of 
the magnetic energy of the intergalactic medium in clusters of galax ies. The source of such 
magn etic fields could be either from the accretion disc of the AGN (jde Kool and Begelman 
I995I ). or f rom dynamo-like effects in the plasmas such as in the solar winds and supern ova 
remnants ( iGiacalone and Jokipii 1 120071 : iKim and Balsara II2006I : iBalsara and Kim 1120051 ). 



We find that magnetic fields of the order of 10 mG, i.e. Alfven velocities of a few 
hundred km s~^, would explain the observed velocity separations between the main and the 
minor components as observed in QSO 2359-1241 and SDSS J0318-0600 (see Bautista et al. 
2010). 
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Another interesting aspect of Eqn. (17) is that it predicts fragments that are faster than 
the main cloud as well as fragments that are slower than the main cloud, as seen in FeLoBALs. 
For example, since W2 = 4: x wi if W2 — Wi = +300 km and v^ave — 1000 km s~^, then 
Vi = 1100 km and V2 = 1400 km s^^. By contrast, if W2 — Wi = —300 km s^^ then 
Vi = 900 km and V2 = 600 km s^^. The faster clouds, in the rest frame of the quasar, 
are expected to outrun the main cloud and remain undisturbed by it. The survival of the 
slower clouds depends on the geometry of the problem. In one dimensional geometry the 
slow clouds would be swept up by the main cloud. In two and three dimensions it is possible 
that magnetic fields produce surface tension in the clouds capable of withstanding the fluid 
bulk pressure. The deceleration of the flow along the axis of travel of the main cloud and 
its fragments yields a bow shock, where the supersonic gas of the main cloud flows smoothly 
around the small fragment (see Shu 1992). 



5. Discussion and Conclusions 

We studied the effects of quasar time variable UV flux on FeLoBAL clouds. In these 
systems absorption from Si II and Fe II comes from a region just before an IF. This IF is 
formed within the absorbing cloud when the cloud crosses a distance Rif from the radiation 
source. At this distance, when the radiation flux drops the ionized region of the cloud will 
retreat at supersonic speeds. A cooling front will follow the ionized region and the IF will 
also retreat and broaden. This creates a pressure imbalance across the IF on time scales 
much shorter than the sound crossing time. Thus, the cooling front yields a rarefaction wave 
that travels into the low temperature region. To conserve momentum across the cloud the 
large pressure imbalance is compensated by an increased speed of the cold gas relative to 
the front. Now, when the flux recovers the neutral rarefled region is re-ionized and re-heated 
by an IF moving away from the ionizing source. This IF can be highly supersonic, driving 
a shock ahead of it and a rarefaction wave in the opposite direction. The shocked gas forms 
new clouds and stabilize at a temperature of ~ 10^ K, like the main cloud. These clouds 
are detached in velocity space from the main cloud and there are rarefaction waves moving 
between the clouds. 

There are two fundamental parameters in determining the speed of the retreating IF. 
One is the ratio of the depth of the cloud to the recombination rate, //r^, which in turn 
is proportional to the column density of the cloud, Nf^. Thus, high column density clouds 
can have faster IFs. The other parameter is the ratio of the variability time scale of the 
ionizing source to the recombination time. Flux variations on scales of ~ 1 year are enough 
for signiflcant effects in FeLoBALs with particle densities of ~ 10^ cm"^, while variabihty 
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on shorter time scales would affect higher density clouds. When the flux drops, broadening 
of the IF slows the cooling front down such that marginally supersonic IF may not lead to 
shocks and waves, but faster IFs will. 

The thermal pressure imbalance across the IFs alone yields velocities to the fragments 
of the order of ~20 km s~^ with respect to the main cloud. Much higher velocities can be 
attained if sizable magnetic pressures are considered. Moreover, in two objects for which 
complete measurements exists, QSO 2359-1241 and SDSS J0318-0600, the observed sepa- 
ration between kinematic components is well reproduced by magnetic field strengths of the 
order of ~ 10 mG. 

Our model predicts that the detached fragments can be either faster or slower than the 
main cloud, as observed in FeLoBAL spectra. Survival of the slower fragments depends on 
the geometry of the problem. We propose that the slower fragments may permeate the main 
cloud if they are clumpy with a small filling factor and/or if the magnetic surface tension of 
the clouds is large enough to make the clouds flow around each other. 

The physical effects presented here may not be exclusive to FeLoBALs, but could have 
implications in a variety of galactic and extragalactic phenomena. For example, in cata- 
clysmic variables, variable stars with circumstellar envelopes such as rj Carinae, and in the 
narrow line regions of AGN. 

We like to thank Nahum Arav and Kirk Korista for fruitful discussions. We also thank 
the anonymous referee for important suggestions. 
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